Relativistic Processes and the Internal Structure 

of Neutron Stars 



D. E. Alvarez-Castillo, S. Kubis 
Institute of Nuclear Physics, Radzikowskiego 152, 31-342 Krakow, Poland 



Abstract. 

Models for the internal composition of Dense Compact Stars are reviewed as well as macro- 
scopic properties derived by observations of relativistic processes. Modeling of pure neutron matter 
Neutron Stars is presented and crust properties are studied by means of a two fluid model. 
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INTRODUCTION 

Neutron stars (NS) are stellar objects that represent one possible end of stellar evolution. 
They are in hydrostatic equilibrium in which the gravitational force is balanced by the 
internal one resulting from the Pauli exclusion principle. Theorized first by Landau 
(1932) and Zwicky (1934), Neutron Stars have been first observed as pulsars in 1967 
by Jocellyn Bell. The latest observations of NS point out that the most extreme physical 
conditions are present in them. They are the densest compact objects (p ~ 10^^ — 10^^ 
g/cm^), fastest rotating stars (as fast as 716 MHz as measured) and fastest moving 
objects in the galaxy (v ~ 1083 km/s). They possess the highest magnetic fields {B — 
10^^ G and largest surface gravity (10^^ cm/s^). Their interior contains superconducting 
material with the highest expected temperature value {Tc = IQP K), and even neutrinos 
can be trapped in proto-neutron stars with temperatures at birth of about 700,000 million 
K[l]. 

NEUTRON STARS MODELS 

Given such extreme conditions various theoretical models have been proposed to de- 
scribe the composition of dense nuclear matter which are translated into an Equation 
of State (EoS). Based on our present knowledge of nuclear interactions, EoSs can be 
categorized in a general form as 

• (Hadronic) Pure neutron matter: mostly n, but also p, e and /i particles. 

• (Hadronic) Hyperonic matter: n, p, hyperons(A, E, S) and leptons (e, /i). 

• (SQM) Strange Quark Matter. Deconfined quarks, where the name goes to the 
heaviest quark that appears as density increases (in this case the strange quark). 

• Boson Condensates {Tt, K): Bosons can be present as a Bose-Einstein condensate 
state of matter (BEC). 



In terms of these EoSs dense compact stars can be classified as: 



- Neutron Stars: both crust and core are described by a hadronic EoS. 

- Hybrid Stars: their crust is hadronic but having a quark matter core. 

- Quark Stars: Only described by a Quark Matter EoS. 

The macroscopic, relativistic structure of a Neutron Star is described as follows. For a 
static, non-rotating star, the Einstein equations give the Tolman-Oppenheimer-Volkoff 
equations [2]: 

dp_ ^ {p+p/c^)G{m + AK?p/c^) dm ^ ^^^^^ 

dr r^{\ —2Gm/rc^) ' dr 

This system requires an EoS of the form p{p) and is to be solved for mass m, pressure p 
and density p inside the star. To find a solution it's necessary to choose a central density 
(Pc) and take into account the boundary conditions that m(r = 0) = 0, m{R) = M and 
p{r = R) = where R and M are the total radius and total mass of the star. As for the 
crust, the moment of inertia / and mass M are derived, in the frame of General Relativity 
[3], as: 

Alcrust - 3 {M^rustR ) i_2GM/7?c2 

where Mcrmt — M —Mcore the difference between the total mass and the mass of the core. 
To determine the latter the crust-core transition point must be known, i.e. the edge of the 
solid crust where the liquid core starts. The EoS of a neutron star is based on nuclear 
models that describe nuclear matter,i.e. a system of interacting nucleons. For the star 
to be stable, charge neutrality (total charge resulting from all its constituents) and beta 
equilibrium (beta and inverse beta reactions taking place at the same rate) must hold. 
Electrons (and at enough high densities muons) form a gas that is negatively charged. 
The crust of the star forms a lattice of nuclear clusters immersed in neutron liquid and 
as such has solid state properties. The core behaves like a liquid and presents a mixture 
of Fermi gases of protons, neutrons and leptons. In the core, the nuclear energy per 
baryon is a function of only baryon number density n and the isospin asymmetry a and 
is defined as follows: 

Enuc (n, a) = Vin) + Sin) * + Q(n) * + ^(a^) (4) 

where a = (n„ — np)/n and n = Up + rin are the neutron and proton number densities. 
Here the last term is negligible since it's contribution is very small. The most interesting 
quantity is the Symmetry Energy (SE) S{n) that has impact in the crust properties of 
the star. It's value at saturation point (the density of the nucleus of an atom), no = 0.16, 
is S(no) = 30 ± 1 MeV, as determined by experiments of isospin diffusion [4] and in 
agreement with the semi-empirical mass formula used to describe nuclei. Values at low 
density have lately been investigated both in theory and experiment [5, 6]. 
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FIGURE 1. Mass vs Radius Relation for EoSs described in Table 1 for non-rotating NS. Blue lines 
represent hadronic EoSs for pure NS. Pink lines are hadronic EoSs that include hyperons. Green lines are 
SQM EoSs. The dashed orange line is the rotational limit imposed by the fastest rotating NS. The dark 
line in the left upper corner represents the causality limit that EoSs must respect. The red band at the 
2 textrmMQ is the measured value for the J1614-2230 NS. More details in [7] from where the figure is 
adapted. 



The core system is then described by thermodynamic quantities derived from this en- 
ergy form (for example p = —dE/dV) at zero Temperature {T = 0) since its contribution 
has no effects in the EoS. To estimate the crust core transition point one may start by 
looking for the instability values against density fluctuations where the system must split 
into two phases i.e. where the compressibility of nuclear matter becomes negative, since 

the condition = > must hold. This marks a lower bound in density values 

where the transition should occur. The line corresponding = is usually called spin- 
odal line. A more elaborated way of addressing this problem is to consider a first order 
phase transition of a two component system. The first consisting only of neutrons while 
the second having neutrons and protons. To ensure stability, mechanical and chemical 
equilibrium must take place by means of the Gibbs conditions: 

J Ji ,,i .,11 ,,i .,11 c;.\ 

P =P , = M« . Me = Me ■ (5) 

Protons are present only in second phase, so for them /i^ > iip. In the two component 
system the average densities and energies are described in terms of the volume fraction 
X occupied by the "I" fluid. They are defined as: 

(e)=ze' + (l-z)e" , {n) =xn' + {\-xW . (6) 

This energy as a function of density is to be compared with the energy of the homoge- 
neous system, the transition point being the value for which the two component system 
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FIGURE 2. Diagram for the two component system for the PALu [8] . Blue contours represent nucleonic 
pressure while orange ones neutron chemical potential. The spinodal line is in red. Green dots are the 
points where the Gibbs conditions are fulfilled. Points to the left of the spinodal belong to the pure neutron 
phase, on the right are points for the nuclear clusters. The crust-core transition takes place at the density 
where the green solid line touches the spinodal. For this model it occurs at ric — 0.0892 fm^^. 

is preferred, where its energy is the lowest. The behaviour of matter close to the crust- 
core transition region is shown in Fig. 2. Any state of matter is represented by a point 
on the X — n plane, where x is the proton fraction x = np/n and n is the baryon number 
density. The green solid line corresponds to the homogeneous neutral matter filling the 
NS core. As the density decreases the homogeneous matter becomes unstable when it 
crosses spinodal line (the red one). All states of matter on left side of the spinodal are not 
stable except the pure neutron matter {x = 0). That means that for low density the system 
splits into two phases: nuclear matter x^O'm the form of finite size clusters immersed 
in pure neutron matter (almost free neutron -i- electron gas). Those two phases are repre- 
sented by two branches of green points for which the coexistence conditions (5) hold. In 
this approach the finite size effect like Coulomb and surface energy are neglected. These 
different effects lead to formation of various structures like rods, plates, and are called 
pasta phases. See [9] for a detailed discussion. 

RELATIVISTIC MEASUREMENTS 

Astronomical observations can shed light on NS macroscopic properties and help to rule 
out theoretical models. The two basic quantities are mass M and radius R. Fig.l shows 



the expected values for families of stars given an EoS. Most processes like Eddington 
flux near the Star's surface, redshift of any luminous signal and spectra from thermal 
bursts involve the compactness parameter ^ where only the ratio of these quantities is 
found [10, 11]. Nevertheless, combining at least two processes is possible to find such 
values, something that has not yet been achieved. For this task, future observations with 
new extraterrestrial telescopes that have been implemented are quite promising. In the 
case of orbiting binary starts only M can be derived with good accuracy in the case of 
orbit reduction (by gravitational wave emission) and Shapiro delay (radar signal delay 
near massive objects) lately measured and resulting in a 2 Mq NS [12]. Macroscopic 
crust properties of NS have also been constrained by pulsar glitching: a sudden spin 
up of the star. For the case of the Vela pulsar a lower bound on the crustal moment of 
Inertia, I/Icmst > 1-4%, has been established by observations based on the superfluity 
model [13]. 

SUIMIMARY 



Neutron Stars are dense compact objects where the most extreme physical conditions 
exist. The macroscopic description being in terms of the General Theory of Relativity. 
For the interior of the star different nuclear models are being used based on terrestrial 
laboratory experiments and extrapolated to the conditions found in NS. Theoretical 
models should be able to reproduce observations, and with the forthcoming data such 
models should start to converge since some of them could be tweaked in their parameters 
to reach the desired values. In particular the latest 2 Mq measurement imposes a strong 
constrain in the EoS. Finally other relativistic processes like the cooling of the proto- 
neutron star could also provide evidence on the interior of NS, since the cooling rate is 
modified by the presence of superfluid matter. All this is reflected in the MvsR relation. 



TABLE 1. 


EoSs used in Fig.l 


Symbol 


References and Specifications 


AkmalPR 


Akmal et al.l998 A18+dv+UIX* (npemu) core, BPS+HP94 outer crust, SLy4 




inner crust. 


SLy4 


SLy4 (npemu) core, BPS+HP94 outer crust, SLy4 inner crust, Douchin and 




Haensel 2001. 


FPS 


BPS below n.drip, then FPS. 


BPAL12 


Prakash 1997. 


BalbNlHl 


BalbNlHl core, SLy4 crust (S.Balberg,October 1997). 


GlendNHS 


Nucleons + Hyperons, Lagrangian mean field theory, Glendenning, N.K. 1985 




ApJ293, 470 . 


SQMa 


SQM EoS: e = p/v^ + St, et = 5.6eo, = 1 /3, Glendenning, August 1990. 


SQMb 


SQM EoS: e = p/v^ + e^, = 5.6eo, = 1, Glendenning, August 1990. 


SQMc 


SQM EoS: e = p/v^ + e^,, e^, = 3.l£b, = 1 /3, Glendenning, August 1990. 
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